Implications of the discrepancy between proton form factor measurements 
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Recent polarization transfer measurements of the proton electromagnetic form factors yield very 
different results from previous Rosenbluth extractions. This inconsistency implies uncertainties in 
our knowledge of the form factors and raises questions about how to best combine data from these 
two techniques. If the discrepancy is due to missing corrections to the cross section data, as has 
been suggested, then the true form factors, related to the proton structure, differ from the form 
factors that parametrize the deviation from point scattering, and different applications will require 
the use of different form factors. We present two extractions of the form factors: a global fit to the 
world's cross section data, and a combined extraction from polarization transfer and cross section 
data. The former provides a parametrization of the elastic electron-proton cross section. The latter 
provides a consistent extraction of the underlying form factors, under the assumption that missing 
terms in the radiative correction explain the difference between the cross section and polarization 
transfer results. 

PACS numbers: 25.30.Bf, 13.40.Gp, 14.20.Dh 



I. INTRODUCTION 

The proton electromagnetic form factors, Ge and Gj\/, 
parametrize deviations from a point particle in elas- 
tic electron-proton scattering, and are related to the 
charge and magnetization distribution of the proton. 
The form factors depend only on Q^, the square of 
the four-momentum transfer, and until recently it was 
believed that the electric and magnetic form factors 
showed approximate scaling, i.e. nearly identical 
dependence H. More recent Jefferson Lab measure- 
ments 0; IS 3 utilized the polarization transfer tech- 
nique to measure the ratio Ge/Gm and found that Ge 
decreases more rapidly than Gm at large Q^. The polar- 
ization transfer measurements are more precise at high 
Q^, and significantly less sensitive to systematic uncer- 
tainties than the Rosenbluth separation measurements. 
However, the two techniques disagree significantly even 
in the region where both yield precise results. 

At the present time, it is not known why the tech- 
niques give different results. The systematic uncertain- 
ties of the polarization transfer measurements, primar- 
ily spin transport and backgrounds, have been carefully 
studied 0. A detailed global analysis of the cross sec- 
tion measurements does not show any inconsistencies 
in the cross section datasets, or yield any likely candidate 
to explain the discrepancy. To resolve the discrepancy, a 
systematic error in the cross section would have to have 
a significant dependence on the virtual photon polariza- 
tion e, = l + 2{l + Q^/4M^)ta.n'^ {Oe/2), where Mp is 
the proton mass and 9e is the electron scattering angle. 
Such a systematic error would have to yield a (5-7)% e 
dependence in the cross section, roughly linear in e, in 
order to resolve the discrepancy. 

There appear to be two possibilities: either a funda- 
mental flaw in the Rosenbluth or polarization transfer 
formalism, or an error in either the cross section or polar- 
ization transfer measurements. Recent works have sug- 



gested that additional radiative correction terms, related 
to two-photon exchange corrections, may lead to an error 
in determining the form factors from the measured cross 
sections If the two-photon exchange mechanism, 

or some other correction that is neglected in the cross sec- 
tion extraction, is the source of the discrepancy, then the 
form factors extracted from a Rosenbluth separation of 
cross section data will not represent the underlying struc- 
ture of the proton, but they will parametrize the elastic 
electron-proton cross section in the usual one-photon ap- 
proximation. Conversely, the true form factors will not 
yield the correct cross sections, and will thus give incor- 
rect results if used as a parameterization of the elastic 
cross section in data analysis. 

If the two-photon exchange term explains the discrep- 
ancy, then the polarization transfer result will relate to 
the true form factors, assuming that the two-photon ex- 
change has a much smaller effect on the polarization 
transfer than on the Rosenbluth extractions. However, 
the existing polarization transfer experiments IS Q 
have extracted the ratio Ge/Gm, rather than the indi- 
vidual form factors. To extract the form factors, these 
data must be combined with cross section measurements 
to determine the absolute magnitudes of Ge and Gm- 
If the two-photon exchange correction modifies the cross 
sections from those calculated from the underlying form 
factors, then it is not possible to consistently combine 
the two kinds of measurements without some assumption 
about the two-photon exchange correction. 

In this paper, we present two extractions of the pro- 
ton form factors. From a global analysis of cross section 
measurements, we extract the "Rosenbluth form factors" . 
From a combined analysis of cross section and polariza- 
tion transfer data, with a "minimal" assumption about 
the nature of the two-photon exchange corrections, we 
extract the "Polarization form factors". If two photon 
corrections are the source of the discrepancy, then the 
Rosenbluth form factors will parametrize the elastic cross 
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section, and are therefore useful as input to analysis or 
simulations that require the electron-proton cross section. 
The Polarization form factors will provide the true form 
factors, which relate to the underlying structure of the 
proton. These form factors are often described as the 
Fourier transformations of the charge and magnetization 
distributions of the proton in the Breit frame, although 
relativistic effects and the fact that each value of cor- 
responds to a different Breit frame lead to substantial 
theoretical difficulties in extracting charge and magneti- 
zation distributions llOl. 



II. GLOBAL ROSENBLUTH ANALYSIS 



The "Rosenbluth form factors" are determined from 
a global fit to elastic electron-proton cross section mea- 
surements. The details of the fitting procedure are de- 
scribed in Ref. [^. For the present analysis we include 
more recent Jefferson Lab measurements of elastic scat- 
tering 0,0,01, as well additional datasets to constrain 
the low behavior 0, 0, 0, 0| to the datasets used 
in Ref. [a, 13- addition, we include all of the high 
Q"^ data, up to 30 GeV^, while the previous analysis was 
limited to 8 GeV^. The older data have updated radia- 
tive corrections, and the small-angle data from Walker et 
al. are excluded, as described in The form factors 
are fit to the following form: 



'P2nQ 
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(1) 

where /ip is the magnetic dipole moment of the proton 
and values are in GeV^. Reasonable fits are achieved 
for > 3. Note that this is a different functional form 
than used in previous fits 0, 0, , which used polyno- 
mials in g = \J Q^. The polynomial in g is a very general 
form, with adequate flexibility to reproduce the data, but 
does not have the proper behavior as 0. 

The fit is quite insensitive to the order of the polyno- 
mial above = 6, except for G^; at large Q^. For 
above 6 GeV^, fits with nearly identical values can 
have Ge/Gm either rise or fall dramatically with Q^. 
This is a result of the reduced sensitivity to and the 
limited e coverage for values above 6 GeV^. To avoid 
unreasonable behavior in the region where G^; is uncon- 
strained by data, we keep the ratio Ge/Gm fixed for all 
values above 6 GeV^. This leads to a fit for Ge which 
is continuous, but not smooth, at = 6 GeV^. Because 
Ge has relatively little contribution to the total cross 
section at these momentum transfers, the cross section 
extracted is still quite smooth, and the value of Gb at 
large values has little effect on the cross section, as 
long as the fit is constrained to avoid \^pGe\ ^ IGaj |. 

The normalization factor for each data set is allowed 
to vary along with the parameters of the fitting functions 
for Ge and Gm- The total from the cross section 
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FIG. 1: (Color online) The "Rosenbluth form factors" (solid 
line) for Ge and Gm relative to the dipole form; Gd = [1 + 
Q'^/Ml]-^, Ml=0.71 GeV^ The dot-dashed line is the pre- 
vious fit to Rosenbluth extracted form factors from |20|l . and 
the dashed curve is the fit to Gjv/fromO, with the form fac- 
tor ratio constrained to give ^pGe/Gm = 1 — 0.13(Q^ — 0.04). 



TABLE I: Fit parameters for the Rosenbluth form factors, 
using the parametrization of Eq. 



Parameter 


Ge (Rosenbluth) 


Gm/^p (Rosenbluth) 


P2 


3.226 


3.19 


Pi 


1.508 


1.355 


Pe 


-0.3773 


0.151 


Ps 


0.611 


-1.14x10"^ 


Pio 


-0.1853 


5.33x10"'' 


Pl2 


1.596x10^2 


-9.00x10"'' 



measurements and normalization factors is: 



E 



Wo: 



(2) 



where <Ti and dui are the cross section and error (exclud- 
ing normalization uncertainties) for each of the N„ data 
points, rjj is the fitted normalization factor for the 
dataset, and drjj is the normalization uncertainty for that 
dataset. We fit to 470 data points (iV<^=443, iVea;p=27) 
with 39 parameters (six parameters each for the electric 
and magnetic form factors, and 27 normalization param- 
eters). 

The result of the global fit to the cross section data is 
shown in Fig.[Tl The fit yields a total of 326.7 for 431 
degrees of freedom, yielding a reduced x^; = l^dof, 
of 0.758. This yields an unreasonably high confidence 
level, indicating that the quoted uncertainties of the mea- 
surements are too large. As was observed in the previous 
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FIG. 2; Gm (top) and Ge (bottom) from direct Rosenbluth 
separation utilizing normalization factors from the global fit. 



fit 1^, the majority of the datasets, 20 out of 27, have 
values of < 1, indicating that most of the experiments 
were overly conservative in estimating their uncertainties. 
Table m lists the parameters for the Rosenbluth form fac- 
tors. The fit includes cross sections for values from 
0.005 to 30 GeV^, and should be valid over this range, 
though the separation of G^; and Gm is only well con- 
strained by the data for <,Q GeV^ . 

The normalization factors were generally smaller 
than the quoted scale uncertainties of the experiments 
(x^ = 18.0 for 27 normalization factors). The average nor- 
malization factor is 0.65%, and the RMS normalization 
factor is 2.7%. The normalization factors are very close 
to those obtained in the previous global fit 0] . The aver- 
age normalization factor differs by approximately 0.5%, 
and the individual normalization factors differ by less 
then 1% for 18 of the 20 experiments. Because the previ- 
ous fit excluded data below (5^=0.6 GeV^, the agreement 
indicates that it is the self-consistency, rather than the 
form factor constraint at = 0, that dominates the 
determination of the normalization factors. 

We can test the self-consistency of the individual data 
sets by comparing the global fit to the results of single- 
experiment extractions of and Gm- By compar- 
ing only the single-experiment extractions, we avoid the 
potentially large and correlated uncertainties that arise 
from the relative normalization of different datasets. 
Comparing the ratio Ge/Gm from the fit to the indi- 
vidual experiments, taken from Refs. [lTlll3l| and the re- 
analysis of older experiments presented in Ref. |H yields 
X^ = 45.3 for 50 data points (x^ = 17.8 for the 20 points 
above Q2 = 1.5 QeV^). 

We can estimate the uncertainties in the form factors 
by performing direct Rosenbluth separations in several 
bins using the full dataset, with normalization factors 



determined from the global fit. For each bin, the data 
are scaled to the average value of the data points in 
that bin, using the global fit as the scaling function. 
bins were chosen so that there are at least three data 
points in the bin, the e range covered is at least 0.3, and 
the correction for scaling each point to the average 
value was <J10% (typically < 2%). The scaling was also 
done using the fits of Refs. and |^], shown in Fig.^ 
Varying the scaling procedure changed the ratios by ^ 
1%, except for the very highest (lowest) points, where 
the change in Ge {Gm) was as much as 3%, but was still 
much smaller than the uncertainty in the extracted form 
factor. 

Figure 13 shows the fits to G^ and Gm, along with the 
direct Rosenbluth separation points, using the normaliza- 
tion factors from the fit. Except for the very low val- 
ues, typical uncertainties on Gm are Ril%, increasing to 
-2% for Q2^10 GeV2 (8% for Q^^SO GeV^). At low Q^, 
the experimental uncertainties become quite large, but 
the constraint on the behavior as ^ yields a much 
smaller uncertainty on the fit. For Ge, the uncertainties 
are (1-2)% at low Q^, but are (5~10)% for intermediate 
values (2-4 GeV^), and grow rapidly as increases. 
Note that the uncertainties in Ge and Gm are highly 
anti-correlated, due to the way the form factors are sepa- 
rated from the cross section measurements. This can be 
seen in the anti-correlation of the deviation of the points 
from the fits in Fig. [3 Thus, the uncertainty on the cross 
sections extracted from this parametrization is not just 
the sum of the uncertainties in the contributions from 
Ge and Gm- Up to Q^k4 GeV^, there is a large body 
of cross section measurements with point-to-point uncer- 
tainties of '^1%. Because the normalization factors are 
determined in the fit, and the residual uncertainty in the 
normalization is small, the absolute cross sections should 
be known to better than 2%. Above (3^=4 GeV^, the 
number of data points decreases, and the uncertainties in 
the cross sections grow, reaching 10% at Q^=25 GeV^. 

Even with the uncertainty related to the discrepancy 
between Rosenbluth and polarization transfer, this fit 
yields a precise parametrization of the elastic cross sec- 
tion in the one-photon exchange formalism. While these 
may not be the underlying form factors of the proton 
(e.g., if there are missing radiative correction terms), this 
is still the appropriate parametrization to use as input to 
a calculation or analysis that requires the elastic cross 
section. Using the form factors derived from the polar- 
ization transfer technique will not yield the correct cross 
section, even in a combined analysis of Rosenbluth and 
polarization transfer such as performed in Refs. 6, 19\. 
More importantly, an inconsistent combination of cross 
section and polarization transfer results can magnify the 
error. Combining a parametrization of Gm from a Rosen- 
bluth analysis with the form factor ratios measured in 
polarization transfer decreases Ge, and thus decreases 
the total cross section, relative to the best fit to the cross 
section data, without allowing a corresponding increase 
in Gm- This leads to form factors which give cross sec- 
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tions that are (4-10)% below the measured cross sections 
at large e over a large range (0.1 < < 15 GeV^). 



III. COMBINED POLARIZATION AND 
ROSENBLUTH ANALYSIS 

While the Rosenbluth form factors yield the best 
parametrization for the cross section in the usual one- 
photon exchange picture, the ratio does not agree with 
the ratio extracted from the polarization transfer tech- 
nique. For the larger values, the polarization transfer 
technique is less sensitive to knowledge of the kinematics, 
radiative correction, and other systematic uncertainties 
that are important in the Rosenbluth separation. 

If this discrepancy is related to a problem in the cross 
section data, then the polarization transfer will yield the 
true ratio of the form factors, but has to be combined 
with cross section data to obtain both Ge and Gm- We 
present in this section a combined analysis of the polar- 
ization transfer and cross section data, which will yield 
the "Polarization form factors" . 

In order to obtain a consistent extraction of the form 
factors, we must make an assumption about the nature 
of the discrepancy. We assume that the difference comes 
from a common systematic error in the cross section mea- 
surements. Analyses of this discrepancy 0, |3| indicate 
that there must be an e-dependent correction of (5-7)%, 
roughly linear in e, for 1 < < 6 GeV^. 

In the combined analysis, we apply a linear, e- 
dependent correction of 6%, to all datasets. This is the 
'minimal' assumption necessary to make the two tech- 
niques consistent, to the extent that a correction that 
was not linear in e, or which modified only some of the 
datasets, would have to be larger. A correction that is 
nearly linear in e and fairly Q^-independent is consistent 
with the form for the two-photon exchange term in the 
analysis of Ref . J} , although the size of the correction in 
Ref. is only ~2%, less than half the size necessary to 
explain the discrepancy. 

We repeat the fit from Sec. ^ but with cross sections 
modified by the linear e dependence, and with the po- 
larization transfer data included in the fit, as described 
in 0. The correction to the cross section could either 
lower the cross section at large e values, or increase it at 
small e values: 

CTci = cro(l - 0.06e) (3) 
CTc2 = fTo(l-0.06(e-l)) = 0.94^,1 (4) 



The first correction is consistent with the form from 
Ref. [g, while the second is consistent with the behav- 
ior of T7|. The second form was chosen for the main 
fit because the correction is small at large e (small (?e), 
where comparisons of positron to electron scattering from 
SLAG [23| set fairly tight limits on the size of two-photon 
exchange. 

The Polarization form factors, from the combined fit 
to the cross section and the 26 polarization transfer data 
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FIG. 3; (Color online) The "Polarization form factors" (solid 
line) for Ge and Gm, relative to the dipole form. The dot- 
dashed line is the previous fit to Rosenbluth extracted form 
factors from 2C|, and dashed curve is the fit to Gm from 
with the form factor ratio constrained to give UpGe/Gm ~ 
1-0.13(Q^ -0.04). 



TABLE 11: Fit parameters for the Polarization form factors, 
using the parametrization of Eq. Q 



Parameter 


Ge (Polarization) 


Gm/hp (Polarization) 


P2 


2.94 


3.00 


Pi 


3.04 


1.39 


Pa 


-2.255 


0.122 


Ps 


2.002 


-8.34x10"^ 


Pio 


-0.5338 


4.25x10"" 


Pl2 


4.875x10"^ 


-7.79x10-*^ 



points from Refs. |^ H, Q is shown in Fig. O The fit 
yields a total of 391.6 for 457 degrees of freedom, xl = 
0.857, including the additional contribution for the 
polarization transfer data (Eq. (8) of Ref. 0). Table Hll 
gives the fit parameters for the Polarization form factors. 

The fit was also performed with the correction of Eq.Ol 
This leads to an overall rescaling of all of the cross sec- 
tions by 6%, relative to the correction of Eq. 0] How- 
ever, this does not yield a simple rescaling of the form 
factors, because each dataset has a normalization fac- 
tor that is determined in the fit, and because the form 
factors are constrained to reproduce the charge and mag- 
netic moment at = 0. While a two-photon correction 
of this size for large e values would appear to be ruled out 
by the SLAG positron-proton measurements 21], an e- 
dependent systematic other than two-photon exchange 
could also resolve the discrepancy. However, this fit 
yields a much worse value: 575.1 for 457 degrees of 
freedom, and so we choose to apply Eq. ^ for our com- 
bined fit. 
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Note that the resuh of the combined fit f Table ITl|l will 
not reproduce the measured elastic cross section in the 
one-photon exchange formalism; it will reproduce the 
modified cross sections of Eq. 0] Therefore, the Polar- 
ization form factors should not be used to model elastic 
electron-proton cross section measurements. However, 
if the 'minimal' assumption that a correction consistent 
with the form of Eq. ^ explains the discrepancy, this 
should yield a consistent extraction of the underlying 
form factors of the proton. 

IV. CONCLUSIONS 

Form factors extracted using the Rosenbluth technique 
provide a parametrization of the deviation of the elastic 
electron-proton cross section from the point-scattering 
cross section. If the cross section has additional correc- 
tions, such as two-photon exchange terms, that are not 
being taken into account, then the Rosenbluth extraction 
does not yield the true proton form factors that relate to 
the structure of the proton. In this case, Ge must be 
extracted from the polarization transfer measurements, 
which yield Ge/Gm, and the cross section data must be 
utilized to determine Gm- While we cannot know how 
to properly combine the polarization transfer and cross 
section data until we understand the cause of the discrep- 
ancy, the uncertainties in Gm that arise from this prob- 
lem are much smaller than those in G^;. The same holds 
true if there is some other correction or combination of 
corrections to the cross section other than the two-photon 
exchange (e.g. Coulomb corrections |2^). It is of course 
possible that the discrepancy is due to a problem with 
the polarization transfer data or technique rather than 
the cross section data. If so, then the Rosenbluth form 
factors represent both the correct cross section and the 
correct nucleon structure. However, there do not appear 
to be any obvious candidates for problems in the tech- 
nique, and the experiment should be less prone to sys- 
tematic uncertainties than the Rosenbluth extractions. 

We have presented two extractions of the proton elec- 
tromagnetic form factors. The Rosenbluth form factors 
come from a global Rosenbluth extraction of the form 
factors from electron-proton elastic scattering measure- 
ments. The Polarization form factors come from a com- 
bined fit to the cross section and polarization transfer 
data, under the assumption that the discrepancy between 
the techniques is caused by a linear, e-dependent cor- 
rection to the cross sections. The Rosenbluth form fac- 
tors give a global parametrization of the elastic electron- 



proton scattering cross section in the one-photon ex- 
change approximation. Even if there is a correction to the 
cross sections, neglected in the one-photon exchange for- 
malism, this parametrization will yield the correct cross 
sections in the one-photon approach. Under the above 
assumption of an unknown correction to the cross sec- 
tions, the Polarization form factors yield the underlying 
form factors, but will not reproduce cross sections, and 
will therefore yield incorrect results if used as input for 
an analysis that requires the elastic cross section, as was 
observed in an analysis of quasielastic scattering from 
nuclei 

Additional data will help shed light on the origin of 
the discrepancy. An i mp roved "Super-Rosenbluth" sep- 
aration measurement [23| completed at Jefferson Lab in 
2002 will yield a precise extraction of Ge/Gm, and will 
determine if the discrepancy can be explained by exper- 
imental problems in the Rosenbluth extractions. A new 
polarization transfer experiment (2jj, approved to run 
at Jefferson Lab, will provide an independent confirma- 
tion of the existing polarization transfer results, as well 
as extending the measurements to higher values. If 
sufhciently improved calculations or direct measurements 
of the two-photon exchange corrections become available, 
we should be able to determine if they are responsible for 
the discrepancy, and if so, remove the current uncertainty 
in combining cross section and polarization transfer mea- 
surements. 

If the discrepancy is explained by two-photon correc- 
tions or some other effect on the cross sections, and we 
have reliable calculations for these effects, then the cross 
section data can be combined with the polarization trans- 
fer data to extract the form factors without ambiguity. 
These form factors will represent the underlying struc- 
ture of the proton and provide a useful parametrization 
of the elastic electron-proton cross section, as long as the 
effect is properly accounted for. Until the discrepancy is 
well understood, however, both sets of form factors are 
necessary, and it is important to use form factors that are 
(1) extracted consistently from the cross section and/or 
polarization transfer data, and (2) appropriate for the 
problem being addressed. 
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